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MIMO TCM WITH CONSTELLATION ROTATION AND 
COORDINATE SWAPPING 

BACKGROUND OF THE INVENTION 
[01] The present invention relates to digital communication, and more particularly to trellis 
coding used with multiple-input multiple-output communication systems. 
[02] Demand for wireless digital communication and data processing systems is on the 
rise. Inherent in most digital communication channels are errors between what was sent and 
what was detected. Such errors are often caused by electrical interference or thermal noise. 
Data transmission error rates depend, in part, on the medium that carries the data; i.e., the 
channel. Typical bit error rates for copper based data transmission systems are in the order of 
one error per 10 6 bits. Optical fibers have typical bit error rates of 1CT 9 or less. Wireless 
transmission systems, on the other hand, may have error rates of 10" 3 or higher. 
[03] The relatively high bit error rates of wireless transmission systems pose certain 
difficulties in encoding and decoding of data transmitted via such systems. A typical coding 
design models errors by an additive white Gaussian noise (AWGN) source added to the 
transmitted signal, partly because of its mathematical tractability and partly because of its 
application to a broad class of physical communication channels. 

[04] One approach to overcoming the effects of noise on data is to encode data at the 
transmitter, in a controlled manner, to include redundancy. The redundancy is subsequently 
used by the receiver to overcome the noise and interference introduced on the signal 
representing the data as the signal is being transmitted through the channel. For example, in a 
simple transmitter, k bits of information are coded by n-bit sequences, where n is greater than 
k, according to some coding scheme. The amount of redundancy introduced by the encoding 
of the data is determined by the ratio n/k, the inverse of which is referred to as the "code 
rate." Codewords representing the n-bit sequences can be generated by an encoder and 
delivered to a modulator that interfaces with the communication channel. The modulator 
might map each received sequence into a symbol. In M-ary signaling, the modulator maps 
each n-bit sequence into one of M= 2 n symbols. Data in other than binary form may be 
encoded, but typically the data is representable by a binary digit sequence. 
[05] To gain bandwidth efficiency and coding, Trellis Coded Modulation (TCM) has been 
developed and combines a multilevel phase modulation signaling set with a trellis coding 



scheme. TCM increases the minimum Euclidean distance between pairs of coded signals to 
reduce the loss from the expansion of signal set and to achieve coding gain with relatively 
simple codes. TCM is described in many places such as in Ungerboeck, "Trellis-Coded 
Modulation with Redundant Signal Sets Part I," IEEE Communications Magazine, Vol. 25, 
February 1987; and in Ungerboeck, "Trellis-Coded Modulation with Redundant Signal Sets 
Part II," IEEE Communications Magazine, Vol. 25, February 1987, both of which are 
incorporated herein by reference for all purposes. 

[06] At the receiving end of a transmission channel, the coded symbols are decoded to 
recover the transmitted data with a goal of minimizing the error rate. The Viterbi algorithm 
is an efficient maximum-likelihood sequence detection method for decoding convolutional 
and trellis coded symbols transmitted over AWGN channels. The Viterbi algorithm is 
described in a number of publications. See for example, Viterbi, "Error Bounds for 
Convolutional Codes and an Asymptotically Optimum Decoding Algorithm", IEEE Trans, on 
Information Theory, Vol. IT-13, April 1967; Forney, Jr., "Maximum-Likelihood Sequence 
Detection in the Presence of Intersymbol Interference" IEEE Trans, on Information Theory, 
Vol. IT-18, pp. 363-378, May 1972; and Forney, Jr., "The Viterbi Algorithm", IEEE 
Proceedings, Vol. 61, pp. 268-278, March 1973, all of which are incorporated herein by 
reference for all purposes. 

[07] In accordance with the Viterbi algorithm, for each received sample, a distance 
between that sample taken at time tj and all the trellis paths entering each state at time tj is 
calculated. In the Viterbi algorithm, the minimum Euclidean distance is selected as the 
optimum branch metric for decoding convolutional and trellis sequences transmitted in 
AWGN channels. 

[08] To increase data rate and operating range without increasing the bandwidth, multiple- 
input multiple-output (MIMO) systems have been developed. In a MIMO system, the 
transmitter includes multiple transmit antennas and the receiver includes multiple receive 
antennas. The transmitter splits the data to be transmitted into a number of streams (typically 
bit streams) or otherwise divides the transmission task and transmits the streams via the 
multiple transmit antennas. The receiver receives the transmitted streams via the multiple 
receive antennas. 

[09] Fig. 1 is a simplified block diagram of a prior art MIMO TCM encoder 10. TCM 
encoder 10 includes a serial-to-parallel converter 12, a convolutional encoder 14, a symbol 
mapper 16 and a serial-to-parallel converter 18. The MIMO system (not shown) of which 
MIMO TCM encoder 10 is a part, includes M t transmit antennas and M r receive antennas (not 
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shown). The channel through which the data encoded by TCM encoder 10 is transmitted is 
characterized by the complex matrix H having the dimensions of M t x M r . The channel is 
assumed to be an AWGN channel. Convolutional encoder 14 in conjunction with symbol 
mapper 14 performs the TCM encoding. 

[10] The labeling convention used herein for alphanumeric symbols represents scalar 
quantities as italic symbols, vectors as lowercase bold symbols, and matrices as uppercase 
bold symbols. Unless otherwise indicated, it is understood that each scalar quantity, vector 
component or matrix element can be a complex number. 

[11] Unencoded data bit streams arrive at serial-to-parallel converter 12 (hereinbelow 
referred to alternatively as converter 12). Out of every u+k bits that serially arrive at 
converter 12, u bits are delivered in parallel to symbol mapper 16 and k bits are delivered in 
parallel to convolutional encoder 14. The k-bits are encoded by convolutional encoder 14, 
which is assumed to have a code rate of k/n. Symbol mapper 16 uses the n encoded bits to 
select one of 2 n cosets of symbols, each having 2 U labels. Symbol mapper 16 then uses the u 
unencoded bits to select one of the 2 U labels from the selected coset. Each label so selected is 
represented by a complex number having a real and an imaginary component thus 
representing a constellation point in the complex plane. Serial-to-parallel converter 18 
converts a serial symbol stream into M t parallel symbols, which are transmitted by the 
corresponding M t transmit antennas. Accordingly, each transmitted vector y has M t complex 
components. Each received vector r has the dimension of M r and is represented by the 
following expression: 

r = Hy + v 

where v is a vector representing noise and/or interference, and is assumed to have a spatial 
auto-covariance matrix, and H represents the channel responses between each transmit 
antenna and each receive antenna. 

[12] Fig. 2 is a simplified block diagram of a prior art MIMO TCM decoder 30 that is 
adapted to decode the data encoded by TCM encoder 10 after this data is transmitted. MIMO 
TCM decoder 30 includes a Logarithmic Maximum A Posteriori (Log-MAP) decoder 30, a 
Viterbi trellis decoder 34, a selector 36 and a parallel-to-serial converter 38. Log-MAP 
decoder 32 is supplied with the received vector r and the channel matrix H and, in response, 
generates an estimate for each of the transmitted symbols; this estimate is also referred to as a 
soft decision of the decoded symbol. In particular, Log-MAP decoder 32 finds the Euclidean 
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distance between the received vector r and the nearest label in each coset for one antenna by 
searching over all other possible transmitted labels on the other M r l antennas and does this 
for each antenna. For each coset j, for all possible labels in coset j and for each transmit 
antenna i, Log-MAP decoder 32 computes a distance metric d(i,j) and a corresponding label 
metric label , as shown below: 

d(i, j) = min (r - HX(£))" A" 1 (r - HX(Ar)) (1) 

k 

label(i, j) = arg min (r - HX(k))" A 1 (r — HX(k)) (2) 

k 

In the above equations (1) and (2): 

H : is the channel estimate matrix, 

A" 1 : is the inverse of the auto-covariance matrix of any of the antennas, 
X : is a M t by CL matrix of the possible transmitted symbol combinations 

on all the M t antennas, where C=2 n is the number of cosets, and L=2 U 

is the number of labels in each coset. 
X(k) : represents all rows of matrix X that have an element of coset k in their 

ith column. 

[13] Metrics d(i 9 j) and label(i,j)so computed are supplied to Viterbi (trellis) decoder 34 
and selector 36, respectively. Viterbi decoder 34 assigns the same distance d(ij) to all edges 
in the trellis that generate an output with elements from coset j. The metrics d(ij) and 
label(ij) corresponding to each transmit antenna i are applied to associated transitions of the 
trellis. For instance, the distance and label metrics associated with transmit antenna 1 are 
applied to the following transitions in the trellis: 

0, M t , 2*M t ... 

Similarly, the distances and label metrics associated with transmit antenna 2 are applied to 
the following transitions in the trellis: 

1, M t +l,2*M t +l,... 
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[14] Viterbi decoder 34 identifies the most likely k unencoded bits received by converter 
12 (see Fig. 1) and the corresponding coset for each channel instance. The most likely coset 
returned from Viterbi decoder 34 is supplied to selector 36. Using the cosets received from 
Viterbi decoder 34, selector 36 identifies the most likely u unencoded bits and the 
corresponding labels for each channel instance. Parallel-to-serial converter 38 receives the bit 
streams associated with the cosets and labels in parallel— as identified by Viterbi decoder 34 
and selector 36— and serializes these bit streams at its output terminals. 
[15] A number of well-known codes that have been devised to increase the reliability in 
data transmission require that errors caused by the channel to be statistically independent. 
This is the case for the AWGN channels. However, there are channels, such as those 
characterized by fading, that exhibit bursty error characteristics. Signal fading often causes 
the signal to fall below the noise level, thus resulting in a large number of errors. One 
conventional method for improving bursty error channels and fading is to interleave the 
codewords in such a way that the bursty channel is transformed into a channel having 
independent errors. In other words, interleaving the codewords before transmission and 
deinterleaving after reception causes burst of channel errors to be spread out in time so as to 
appear as random errors. 

[16] It is desired to have a technique for increasing the reliability with which data is 
transmitted through fading channels in a MIMO TCM system. 



BRIEF SUMMARY OF THE INVENTION 
[17] In accordance with the present invention, in a multi-input multi-output 
communication system having M t transmit antennas and M r receive antennas, complex 
symbols are arranged so that distinct symbols have distinct real and imaginary components 
for some set of symbols, such as by being rotated by an angle. The real components of each 
associated pair of arranged symbols are then used to form an input symbol stream- 
representing a complex number— that is subsequently transmitted on, e.g., the first one of the 
transmit antennas. The imaginary components of each associated pair of arranged symbols 
are used to form a second input symbol stream— representing another complex number— that is 
subsequently transmitted on, e.g., the second one of the transmit antennas. Accordingly, one 
of the real and imaginary components of each symbol is swapped at the transmitting end. In 
some embodiments, each arranged symbol is interleaved before being swapped. 
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[18] If the number of transmit antennas is equal to two, (i.e., Mt=2) the input bit stream x x 
associated with encoded symbols Encj and Enc2 and transmitted on, e.g., antenna 1 is as 
follows: 

x , = Re(Enc { ) + j * Re(Enc 2 ) 

Therefore, antenna 1 transmits a complex number that is formed by the real components of 
encoded symbols Enci and Enc 2 . Similarly, the input bit stream x 2 associated with encoded 
symbols Enci and Enc 2 and transmitted on, e.g., antenna 2 is as follows: 

x 2 = Im(Enc 1 ) + j*Im(Enc 2 ) 

Therefore, antenna 2 transmits a complex number that is formed by the imaginary 
components of encoded symbols Enci and Enc 2 . 

[19] At the receiving end, a decoder receives a vector r of the transmitted bit streams 
—received on the M r receive antennas— and a channel matrix H that includes M t rows and M r 
columns. For each transmit antenna i (i is an integer varying from 1 to M t ), for all L possible 
labels in each coset j, a Log-MAP decoder at the receiving end computes an associated 
distance d(i,j) as shown below: 

d(i 9 j) = min(r - h,x, (*) - H^.x)" A" 1 (r - h {Xj (*) - H W9t/ x). 

where xj(k) is the k th constellation point in coset j for antenna i, h x is the i th column of H, 
H n;ti are the remaining columns of H, A" 1 is the inverse of the antenna auto-covariance 
matrix, X is a (M r l) by CL matrix of the p==2 (u+nr(Mt ' 1} possible transmitted symbol 
combinations on the remaining antennas, C=2 n is the number of cosets, and L=2 U is the 
number of labels in each coset. 

[20] Next, using the channel estimates H and r, for each transmit antenna i, for all L 
possible labels in each coset j, the Log-MAP decoder computes an associated label 
label(i,j) as shown below: 
[21] 
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label{U j) = arg min (r - h,x , (*) - U^xf A _1 (r - h,x y (*) - H fl ,x) 

[22] A Viterbi decoder receives the distances d(i,j) and labels label(ij), and in response, 
identifies the most likely transmitted coset. Subsequently, a selector selects the transmitted 
label belonging to coset identified by the Viterbi decoder. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[23] Fig. 1 is a simplified block diagram of a MIMO TCM encoder, as known in the prior 
art. 

[24] Fig. 2 is a simplified block diagram of a MIMO TCM decoder, as known in the prior 
art. 

[25] Fig. 3 is a simplified block diagram of a MIMO TCM encoder, in accordance with 
one embodiment of the present invention. 

[26] Fig. 4 is a simplified block diagram of a MIMO TCM decoder, in accordance with 
one embodiment of the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 
[27] Fig. 3 is a simplified block diagram of a MIMO TCM encoder 100, in accordance 
with one embodiment of the present invention. MIMO TCM encoder 100 includes, in part, a 
serial-to-parallel converter 102, a convolutional encoder 104, a symbol mapper 106, a serial- 
to-parallel converter 108, and a coordinate swapper 110. Data supplied by coordinate 
swapper 1 10 is transmitted by at least two transmit antennas (not shown) and received by M r 
receive antennas (not shown). The channel through which the data encoded by TCM encoder 
100 is transmitted is characterized by the complex matrix H having the dimensions of M t x 
M r . In the following description of MIMO TCM encoder 100, it is assumed that M t is equal to 
two. It is understood, however, that the present invention applies to MIMO TCM encoders 
having more than two transmit antennas. 

[28] Unencoded data bit streams arrive at serial-to-parallel converter 102 (hereinbelow 
referred to alternatively as converter 12). Out of every u+k bits that serially arrive at 
converter 102, u bits are delivered in parallel to symbol mapper 106 and k bits are delivered 
in parallel to convolutional encoder 104. The bits do not need to be delivered exactly in 
parallel, but can be treated logically as parallel bits. The k bits are encoded by convolutional 
encoder 104 which has a code rate of k/n. Symbol mapper 106 uses the n encoded bits to 
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select one of 2 n cosets of symbols, each having 2 U labels. Symbol mapper 106 then uses the u 
unencoded bits to select one of the 2 U labels from the selected coset. Each label so selected is 
a complex number having a real and an imaginary component. The serial bit streams 
associated with the real and imaginary component of each selected symbol are then rotated by 
symbol mapper 106 so that the real and imaginary components of the resulting symbols have 
distinct coordinates; standard constellation rotation techniques may be employed. 
Alternatively, any other technique for making the real and/or imaginary components more 
distinguishable might be used. 

[29] Constellation rotation (for clarity of description, the distinguishing of components is 
assumed to be constellation rotation in this portion of the description) allows every symbol in 
the constellation to be uniquely identified based on either its real or imaginary component. 
The rotated serial bit streams are then coordinate interleaved by symbol mapper 106. 
Coordinate interleaving results in the real and imaginary components of a symbol being 
transmitted by different antennas. Used in conjunction with constellation rotation, coordinate 
interleaving increases the diversity order of the system by providing additional redundancy. 
Specifically, if the signal from one transmitter is not received (due to channel fading), only 
one of the real and imaginary coordinates of a symbol is lost. The symbol may still be 
determined based on the other coordinate, which is received from a different transmitter that 
has different fading properties. 

[30] Serial-to-parallel converter 108 converts the rotated and interleaved serial bit streams 
associated with each selected symbol to M t parallel bit streams each having a real and an 
imaginary component. Each of the M t antennas (M t =2), in accordance with the present 
invention, transmits a signal representing a complex number formed by either only the real 
components of the selected symbols or the imaginary components of the selected symbols. 
For example, assume that encoded symbols Enci and Enc 2 are selected for transmission. Each 
of the encoded symbols Enci and Enc 2 includes a real component and an imaginary 
component. In accordance with the present invention, the input bit stream to, e.g., antenna 1 
is as follows: 

x , = Re(Enc 1 ) + j * Re(Enc 2 ) 
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[31] In other words, antenna 1 transmits a complex number that is formed by the real 
components of encoded symbols Enci and Enc 2 . Similarly, in accordance with the present 
invention, the input bit stream to, e.g., antenna 2 is as follows: 

x 2 = Im(Enc, ) + j * Im(Enc 2 ) 

Therefore, antenna 2 transmits a complex number that is formed by the imaginary 
components of encoded symbols Enci and Enc2. 

[32] Fig. 4 is a simplified block diagram of a MEMO TCM decoder 200, in accordance 
with one embodiment of the present invention. MIMO TCM decoder 200 includes, in part, a 
Log-MAP decoder 202, a Viterbi decoder 204, a selector 206 and a parallel-to-serial 
converter 208. MIMO TCM decoder 200 is adapted to decode symbols encoded via an 
encoder, such as MIMO TCM encoder 100 shown in Fig. 3, and transmitted via a channel 
with a characteristic estimate matrix H. 

[33] In the following, it is assumed that encoded, rotated, interleaved and swapped symbol 
pair components forming associated bit streams are transmitted on M t =2 transmit antennas 
and are received on M r receive antennas. Therefore, each transmitted symbol vector y has M t 
complex components. It is further assumed that each received vector symbol r has the 
dimension of M r and is represented by the following expression: 

r = Hy + v 

where v is the noise and/or interference vector which is assumed to have a spatial auto- 
covariance matrix. 

[34] Log-MAP decoder 202 receives the channel matrix H and vector r, and in response, 
supplies to Viterbi decoder 204 and selector 206 a distance metric (alternatively referred to 
hereinbelow as a distance) and a label metric (alternatively referred to hereinbelow as a label) 
associated with each transmitted symbol. For each transmit antenna, a set of associated 
distances and labels are computed for each m th channel instance given the channel estimates 
H(m) (alternatively referred to hereinbelow as H) and the received vector r(m) (alternatively 
referred to hereinbelow as r). A channel instance is associated with one of the M t transmit 
antennas and one of the M r receive antennas, so the m th channel instance is one of M t x M r 
channel instances. 
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[35] The distance between the received vector r and the nearest label in each coset is 
determined for each antenna by searching over all other possible transmit symbols sets on the 
other M r l antennas. In other words, for each transmit antenna i (i is an integer varying from 
1 to M t ), for all L possible labels in each coset j, a distance d(i 9 j) is calculated as shown 
below: 



d{U j) = min(r -h §Xj (*) -H nH x) H A" 1 (r - h,.* , (*) - H^x). 

where: 

Xj(k)\ is the k th point in coset j for antenna i 

h; : is the i th column of H, 

H n?ii : are the remaining columns of H, 

A" 1 : is the inverse of the antenna auto-co variance matrix, 
X : is a (M t - 1 ) by CL matrix of the p =2 (u+n) *< Mt " 1 > possible 

transmitted symbol combinations on the M r l other antennas, 
where C=2 n is the number of cosets, and L=2 U is the number of 
labels in each coset. In other words, X includes all the possible 
transmitted symbols except those transmitted on the i th antenna. 

[36] Each xj(k) and each component of matrix X is rotated, coordinate interleaved and 
swapped as described above. Accordingly, Log-MAP decoder 202 is adapted to deinterleave, 
unrotate and unswap the data streams associated with each received symbol. Edges in the 
trellis that generate an output with elements from coset j, can be assigned the same distance 
d(i,j) . The distances are computed for each of the other transmit antennas, i=2,. . .M t . Next, 
using the channel estimates H and r, the corresponding labels are computed for each edge in 
the trellis decoder, as shown below: 

label{U J) = arg min (r - h,*,(*) - H^xf A" 1 (r - h iXj (*) - H_x) 

k 

[37] The distances d(ij) and labels label(ij) corresponding to each transmit antenna are 
applied to associated transitions of the trellis. For instance, the distances and labels associated 
with the first transmit antenna (i.e., i=l) are applied to the following trellis transitions: 
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0, M t , 2*M t . 



Similarly, the distances and labels associated with the second transmit antenna (i.e., i=2) are 
applied to the following trellis transitions: 



l,M t +l,2*M t +l,... 



[38] Viterbi decoder 204 receives the distances d(ij) and labels label(ij) generated by 
Log-MAP 202. In response, Viterbi decoder 204 identifies the most likely transmitted coset. 
Using the cosets identified by Viterbi decoder 204 and the labels generated by search- 
efficient LOG-MAP 202, selector 204 selects the transmitted label belonging to coset 
identified by Viterbi decoder 204. Parallel-to-serial converter 208 receives the cosets and 
labels supplied by Viterbi decoder 204 and selector 206 in parallel. In response, parallel-to- 
serial converter 208 supplies these streams serially at its output terminals. 
[39] It is understood that the various blocks in encoder 100 and decoder 200, in accordance 
with the present invention, may be implemented by software code executed by a central 
processing unit, by dedicated hardware or by a combination of hardware and software. 
[40] The above embodiments of the present invention are illustrative and not limitative. 
The invention is not limited by the number of deployed transmit or receive antennas. The 
invention is not limited by the number of cosets or the number of labels within each coset. 
The invention is not limited by the number of bits used to select a coset or the number of bits 
used to select a label from a coset. Other additions, subtractions or modifications are obvious 
in view of the present invention and are intended to fall within the scope of the appended 
claims. 
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